Phase transitions in topological lattice models via topological symmetry breaking 
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We study transitions between phases of matter with topological order. By studying these tran- 
sitions in exactly solvable lattice models we show how universality classes may be identified and 
critical properties described. As a familiar example to elucidate our results concretely, we describe 
in detail a transition between a fully gapped achiral 2D p-wave superconductor (p + ip for pseu- 
dospin up/p — ip for pseudospin down) to an s-wave superconductor which we show to be in the 2D 
transverse field Ising universality class. 



Motivated in part by growing interest in topological 
quantum computation, a considerable effort has been in- 
vested in understanding systems which realize topologi- 
cally ordered phases of matter. Though much is under- 
stood about the topological properties of these phases 
and their possible applications to quantum computing [I]- 
httle is presently known about transitions between 
proximate phases of different topological order. While a 
powerful formalismf^", known as "topological symmetry 
breaking" (TSB), has been worked out to identify when 
pairs of phases can in principle be connected by a con- 
tinuous transition, the critical properties of these tran- 
sitions have not been systematically studied. Recently 
increasing interest has been focused on such phase transi- 
tions, highlighting both their importance to understand- 
ing physical systems such as quantum Hall bilayers[5], 
and their potential as theoretical models displaying crit- 
ical behavior unlike that of more conventional statistical 
mechanical systems [6 . 

In this paper we take a step towards building a sys- 
tematic understanding of the relationship between TSB 
transitions and second-order phase transitions involving 
a broken symmetry. Our main result is to describe a 
large class of transitions between distinct, topologically 
non-trivial phases, and show that the effective theory of 
these transitions is exactly that of the 2D transverse field 
Ising model (TFIM). Specifically, we study such phase 
transitions in the lattice models of topological matter in- 
troduced by Levin and Wen[7j- We show that for models 
described by an SU{2)k x S'J7(2)j, topological field the- 
ory, we can add a perturbation to the Hamiltonian which 
condenses a bosonic field, and that at energy scales be- 
low the quasi-particle gap we can map the perturbed sys- 
tem exactly to a transverse field Ising model. This map- 
ping reveals the precise relationship between the phase 
transition in the topological system and a related global 
symmetry-breaking transition. The lattice models we 
construct also give interesting proof-of-principle exam- 
ples of the TSB scenario pi, in which we can track explic- 
itly the fate of the topological order after condensation. 

While the lattice models we study are best viewed as 
toy models, the long-wavelength behavior of these toy 



models should be applicable in analogous real physical 
systems. Consequently we will center our discussion 
around the example probably most familiar from the lit- 
erature, where the topological order is that of a chiral 
p-wave superconductor [5]. Specifically, we study a Levin- 
Wen model describing a bi-layer net achiral p-wave su- 
perconductor, in which the order parameter is p -I- for 
pseudospin- up (top layer), and p—ip for pseudospin down 
(bottom layer). The Ising phase transition we describe 
in this case is to a phase which is topologically an s-wave 
superconductor. 

Let us first describe briefly the exactly solvable Levin- 
Wen lattice Hamiltonians for the honeycomb lattice[7]. 
For simplicity, we focus on the lattice model realizing 
an SU{2)2 x SU{2)2 topological order, which is that of 
the achiral p + ip/p — ip superconductor. In this case, 
each edge in the lattice model can be in one of three 
states, which we label 1, cr, and ip. Two types of operators 
act on these states: string operators s,-, which raise and 
lower all edge labels along the trajectory of the string s, 
and vertex projectors By- In both cases, their action is 
determined by the fusion rules of SU{2)2 Chern-Simons 
theory (^xip — 1, ax^ — a, and a x a ~ 1 + ip) which 
describes the topological excitations in the chiral p-wave 
superconductor. By projects onto states in which the 
three labels entering the vertex V appear in one of the 
following 'allowed' combinations: 




(1,1,1) {ir,cr,4i) (o-,(i,l) 



which are chosen such that the three labels incident at 
each vertex can fuse to give the identity in SU{2)2- The 
string operators act on edge labels according to 

s„\l)^\a) sM)^\a) sM^\l) + \iP) 

s^\<y)^\a) s^m-^\l) (2) 

Here the coefficients of proportionality are dictated by 
properties of SU{2)2 theory, as described by Ref. [71 but 
will not be crucial for our discussion. 
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The Hamiltonian is 

i? = -^(l-Sy)-^(l-<) (3) 
V P 

where ~ \{^~^ V^s^ + s^) is a projector com- 

posed of string operators acting on the edges border- 
ing the plaquette P. Since aU of the projectors in ([3| 
comniute[7j, the spectrum of H can be determined ex- 
actly. Excitations consist of 8 types of quasi-particles 
CTL.i?, ^L,R, ctlctr, ctliI^r, ipLCTR, and ^^^L■^^^R. The sub- 
scripts R and L denote the chirahties of the particles. 

We may understand the spectrum of this system by 
considering the corresponding excitations in the super- 
conducting system. The p + ip/p — ip superconductor 
should be thought of as two independent layers with op- 
posite (L, R) chirahties. <Jr^l corresponds to the super- 
conducting vortex (together with its associated 0-energy 
Majorana fermion bound state) in the top (i?) or bot- 
tom (L) layer; ^r,l are fermionic quasiparticles in the 
two layers, and the other 6 excitations are composites 
of fermions and/or vortices in both layers. The mutual 
statistics of particles in the same layer are those of the 
chiral p-wave superconductor whereas excitations in dif- 
ferent layers have trivial mutual statistics. Note also that 
the tjj particles have Z2 symmetry, meaning they fuse 
with themselves to give 1. In the superconducting pic- 
ture this reflects the fact that fermions are conserved only 
mod 2. 

In this model, the particle = ipL^^R is a boson and 
can therefore in principle condense in a second-order 
TSB-type phase transition. Here we study this transi- 
tion, showing that as i/i is a Z2 field, the transition is in 
the TFIM class. 

To condense (j) i^i the lattice model, we modify the 
Hamiltonian (|3| in two ways. First, noting that (j) vi- 
olates only plaquette terms [9], we decrease the gap to 
creating (/> by modifying the plaquette projector: 

= 1 (4°) + B(f )) + (4°) - 4«) (4) 

where i?^"* — | (l — V^Sg- + s^) is a projector with 
eigenvalue 1 if the plaquette P contains a particle, and 
otherwise. Second, we add a term 

Hi = -j.Y.^-ir'^'^ (5) 

e 

where ncr(e) is 1 if edge e carries label a and is oth- 
erwise. It is easy to check that the operator (— 
applied to the ground state creates (j) particles on the 
pair of plaquettes bordering e (or hops an existing (p from 
one plaquette to another). Hi is equivalent to adding an 
energy penalty to all edges carrying the a label. In the 
superconducting analogue, the creation of a = ipLfpR 
on two neighboring plaquettes is quite a natural term to 
add to the Hamiltonian - corresponding to a four-fermion 



interaction term that creates or annihilates cooper pairs 
in the pseudospin (inter-layer)singlet channel. 

To describe the resulting phase transition, we exploit 
the fact that both Hi and B^p^ commute with By at 
every vertex, and neither creates plaquette violations 
other than (j). We can then study the phase transi- 
tion in the reduced Hilbert space consisting only of the 
ground state plus some number of particles created. 
Since (/) is a Z2 field, we describe the states in the low- 
energy sector of the lattice model by an Ising variable 
ri0 = i {Sz -1-1) = 0, 1 on each plaquette (together with 
a label a identifying the ground state sector if the system 
is on the torus). The vertex projectors, as well as the pla- 
quette term ^ (^Bp^ + i?^"*^, act as the identity on this 
space. Within this low-energy sector, the plaquette term 
4f (^Bp^ ~ B^p^^ acts like the spin operator JzS'z, diago- 
nal in the basis, while Hi acts like — X]<ij> Sx^ Sx'' 
as it flips pairs of spins on neighboring plaquettes. One 
can check that the operators in Eqs. Q [s]) satisfy the 
appropriate commutation relations. 

We thus arrive at an effective theory for the phase tran- 
sition which is exactly the transverse field Ising model, 
with spins on sites of the dual lattice. The paramag- 
netic phase, in which dominates, corresponds to the 
initial SU{2)k x S'J7(2)^ topological phase; the ferromag- 
netic phase, in which Jx dominates, contains an indefinite 
number of (j) particles in each plaquette, and corresponds 
to the condensed phase. The transition between these 
two phases is simply the quantum phase transition of the 
2D TFIM. 

At this juncture, the reader may wonder how it is that 
we have mapped a transition which changes the topologi- 
cal order onto one which breaks a global symmetry. Since 
topological order is associated with long-ranged statisti- 
cal interactions (or generalized Berry's phases) which can 
only occur in the presence of gauge fields, this is indeed 
a surprising result. One way to understand it is to note 
that is essentially the vortex of a Z2 gauge theory; the 
special form of the Hamiltonian we use here ensures that 
no objects behaving like electric sources of this gauge field 
are present in the ground state anywhere in the phase 
diagram. It is because of this that we may exploit the 
duality between the TFIM and Ising gauge theory[TU] to 
study the phase transition. For a more generic choice 
of Hamiltonian this mapping will no longer hold, as the 
appropriate dual theory would again be an Ising gauge 
theory. However, it is known matter sources do not affect 
criticality in the Ising gauge theory [llj - hence we expect 
that the transition will still be described by an Ising crit- 
ical point away from the solvable limit considered here. 

Having identified a mechanism for condensing the par- 
ticle (f), we now describe the topological order of the con- 
densed phase. We will see that the condensed phase of 
the lattice model can be mapped exactly onto Kitaev's 
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Toric code (TC) [5] , which is known to be the same topo- 
logical order as that of an s-wave superconductor |12j . 

To understand the topological order of the condensed 
phase, it is helpful to consider the limit — 0, in which 
the lattice Hamiltonian is again exactly solvable. Here 
the plaquette projector takes on a particularly simple 
form: Bp — 1 + s^, which leaves ria- on each edge un- 
changed. When Jx > 0, this implies that the exact 
ground state contains no edges labeled a, as Hi induces 
an energy cost 2Jj. for each of these. In this limit, the 
label a is confined and we drop it from the low-energy 
theory. This leaves a solvable lattice model containing 
only the two edge labels 1 and ■0, with the Hamiltonian 
of the Toric code (TC)j2] on the honeycomb lattice. 

Our next task is to understand how the 8 excitations 
of the initial model are related to the TC's 3. To identify 
the spectrum of the condensed phase, we must account 
for three effects |4|. First, ipL and ipp are mixed by fu- 
sion with the condensate, and become indistinguishable. 
Such identification of particles is a generic feature of con- 
densation transitions. In the language of the supercon- 
ductor, once a (c^c^) becomes nonzero, fermion number 
is not conserved mod 2 in each layer independently, but 
only in the overall system as a whole (i.e., Z2 x Z2 is 
broken to a single Z2). Second, the label a becomes 
confined, effectively eliminating the four quasi-particles 
(JLt<^r,'^l<^Ri and ctlV'-R from the spectrum. In the su- 
perconducting picture, this is because a vortex in a single 
layer engenders an energetically costly 'branch cut' in the 
s-wave condensate. Finally, TSB predicts that the bound 
state (Jl'^r splits into two distinct types of excitation in 
the condensed phase|4j. We will show that in the lattice 
model this splitting can be understood as a result of a 
new conservation law. 

The energetics of confinement are easily understood 
in terms of the Ising model representation. The opera- 
tor , which measures the flux of (j) through a plaque- 
tte P, is simply \ {sf^ - sj?^) = y/2sa{P) ~ the cr- 
labeled raising operator applied to all edges surrounding 
this plaquette. In the absence of vertex violations, the 
cr label always forms closed loops on the lattice, which 
constitute domain walls in the S*^ basis, since we identify 
(_xyia(u) ^ S^{i)S'^{j). In the ferromagnetic phase of 
the Ising model, these domain walls are confined. The 
energy cost of a string of edges labeled a is thus linear 
in the string length, and the label a disappears from the 
long-wavelength theory. 

Alternatively, we may exploit the duality between the 
TFIM and Ising gauge theory, given by 

i borders p 

(6) 

where i is the edge between the two plaquettes i and j. 
This maps the domain walls of the Ising model - given 
in our lattice model by edges where n„ = 1 - onto the 



electric field = — 1 of the Ising gauge theory. The 
symmetry-breaking transition in the TFIM is thus a con- 
fining transition in the gauge theory [13]. This allows us 
to conclude that open a strings (which join pairs of a^ ji 
vortex excitations) are also confined. Readers might ob- 
ject that this identification of a strings with the Ising 
electric flux is not exact, as ul^r are not bosons. This 
disagreement in statistics is, however, immaterial to the 
energetics of confinement, for which we require only static 
sources. 
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FIG. 1. The action of the operator So-cr which creates the 
quasi-particle OL<yR which splits in the condensed phase. The 
product So-ji can be resolved in terms of strings X,Y = 
which act only on the edge labels, and a phase (depicted 
diagramatically by the ring labeled cr) each time the string 
crosses an edge in the lattice. (On edges without crossings, 
the phase must be 1). In the uncondensed phase X,Y = 1 
ijj on each edge; after condensation X = Y and there are two 
distinct excitations X = Y = ^l) and X = Y = 1. 

To trace the origins of the splitting, we first must un- 
derstand the form of the operators which create aL<JR- In 
the p + ip/p — ip superconductor, this is a pair of vortices, 
one in each layer, which are bound together. This exci- 
tation may appear in one of two states: the total fermion 
number of the pair of vortices may be 1 or ip. However, 
this fermion number will change if a single ctl h vortex in 
one of the layers braids with this bound pair, and is not 
a topologically protected quantity. Hence it is conserved 
in the phase with s-wave pairing, where aL^n vortices are 
confined, but not in the pure p-wave phase. 

In the lattice model, the same effect is seen in the 
quasi-particle creation operators. The operator Saj^s^^^ 
which creates a^afi particles acts with a combination of 
a phase factor at each edge crossed by the string opera- 
tor, and by acting with si_^ G {si, s^} on each edge that 
s runs along. (This is depicted diagramatically in Fig- 
ure [T]). The label 1,^/' carried by this string corresponds 
to the net fermion number of the vortex pair discussed 
in the previous paragraph. Table |l] gives the action of 
the components of Sa-a- Here again we find that So-o- oc- 
curs in two distinct flavors, sj.^ and sf^, distinguished 
by whether or not they interchange the edge labels 1 and 
-0. Before condensation only the symmetric combination 
of these comprises a topological excitation (other com- 
binations have an energy cost which grows with string 
length). After conflning the a strings, however, we flnd 
two distinct quasi-particles, aai and aa^. 

To summarize, in the condensed phase we flnd three 
quasi-particles ip^aai, and cra^. These correspond ex- 
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TABLE I. The action of So-o- on an edge label i has two 
components: the action of a string sx which maps i to i' , and 
a phase 6. (Combinations not shown give 0.) X and Y may 
be distinct only if the edge i crossed by the string operator 
carries the label a. After condensation a is confined, and the 
value of X is conserved along the string. 



actly to the three excitations in the TC: the two bosons 
aai and aa^ are the vortex m and electric source e re- 
spectively; the fermion ip is the combination em. The 
phases inherited from the original ctlO'b. creation oper- 
ators ensure that braiding any two different particles 
around each other incurs a phase of —1. 

To clarify the fate of the topological order across the 
phase boundary, we consider the ground state degeneracy 
on the torus. In the SU{2)2 x SU {2)2 phase there are 9 
ground states \^a), which can be identified by the flux a 
through one of the non-contractible curves on the torus. 
These fall into two classes, distinguished by the opera- 
tors Li = He on Ci(~l)""^''^ where Ci,i = 1,2 run along 
edges in the dual lattice around the two non-contractible 
curves on the torus. Four of the ground states have at 
least one Lj = —1. Because a can only appear in closed 
loops, this means that these states must have a string of 
edges labeled a which run across the width of the system 
in one of the two directions; in the confined phase, where 
such a string incurs an energetic cost linear in the system 
size, they are no longer ground states. Of the remaining 



five, the two anti-symmetric combinations (|f2 



and i^xji — vanish in the condensed ground 

state. Hence these also become gapped and split off from 
the ground-state sector, leaving only the two symmet- 
ric combinations + I^V'j?) ^'^d (I^^Vi.'/'h) + l^i))- 
Finally, as ip appears only in closed loops after condensa- 
tion, jrio-iCTj;) splits into two states, distinguished by the 
eigenvalue of He on ^.^ (— 1)"'''^'^^. This gives four ground 
states on the torus, as expected for the TC. 

The transition studied so far, between SU {2)2 x SU{2)^ 
and the TC, is prototypical of a class of transitions be- 
tween pairs of topological phases, all of which exhibit 
transverse-field Ising criticality. We postpone a detailed 
description of these other transitions to a future work; 
here we merely outline their general features: Consider 
a Levin- Wen Hamiltonian describing the topological the- 
ory SU{2)k X SU{2)i^, which has quasi-particles of spin 
0, |, in two mutually non- interacting sectors R and 
L. We condense the boson (pk = 1^1^ by modifying 
the plaquette projector as in Eq. Q, and adding a bo- 
son creation term where s(e) is the spin on 
edge e. As (pk is always a Z2 field which commutes with 
the vertex projectors, the dynamics of the phase transi- 



tion is described by the transverse field Ising model. All 
half-odd integer spin edge labels map to the domain wall; 
hence particles of net half-odd integer spin (of the form 
^L^^2^fl' ^^""^ versa) are confined in the condensed 
phase. As before, pairs of particles (*Ljfl, ^-j^^^-j^^ 

which mix by fusion with ^ are no longer distinct in 

,2 I, 

the condensed phase. This leaves ^ + | + 1 potential 
quasi-particles. 

As in the k — 2 case, it remains to ask whether any 
of these quasi-particles will split. If k is odd, there is no 
splitting. This leads to a theory whose topological struc- 
ture is SO{3)k X 50(3)^. - meaning that the only effect of 
condensation is to eliminate half-odd integer spin excita- 
tions from the theory. The final theory is simple in this 
case because the initial particle spectrum is of the form 
{ai,a2, ■■■ar/2',4' ^ ai,(p x a2,---4> x 0^/2}, and conden- 
sation merely confines particles of net half-odd integer 
spin, and identifies pairwise excitations ^'^2"^ j;, ^^2"^ ^ = 
^'"^i"^ i^ ^^^i^^ ji- If k is even, the criteria of Ref. H 
dictate that the particle splits into two compo- 

nents. In the lattice model, these are distinguished by 
whether the associated edge string operator carries even 
or odd integer spin. In this case the final theory is not 
a doubled theory with decoupled right- and left- handed 
sectors, but a more involved type of achiral theory known 
as a Drinfeld double P3]- 

In fact, the protocol outlined here can be extended to 
a wider class of transitions for models where an achiral 
boson is condensed, resulting in a transition of the 
q-state transverse field Potts type. 

Here we have described a phase transition between 
phases exhibiting the topological orders oi &p + ip/p — ip 
superconductor, and that of an s-wave superconductor 
(or equivalently, the Toric codelH [H]). By realizing 
the initial phase in a solvable Levin- Wen [7j type lattice 
model, where the phase transition can be induced by a 
simple deformation of the Hamiltonian, we show that and 
that the phase transition is that of the 2D transverse-field 
Ising model. The analogous transitions for the topologi- 
cal order of SU{2)k x SU{2)^, can also be studied in the 
lattice model, leading to a family of condensation transi- 
tions which are all in the same Ising universality class. 
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